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ABSTRACT 

The upcoming new generation of spectroscopic galaxy redshift surveys will provide large samples of cosmic 
voids, the distinct, large underdense structures in the universe. Combining these with future galaxy imaging 
surveys, we study the prospects of probing the underlying matter distribution in and around cosmic voids via 
the weak gravitational lensing effects of stacked voids, utilizing both shear and magnification information. The 
statistical precision is greatly improved by stacking together a large number of voids along different lines of 
sight, even when taking into account the impact of inherent miscentering and projection effects. We show that 
Dark Energy Task Force (DETF) Stage IV surveys, such as the Euclid satellite and the Large Synoptic Survey 
Telescope (LSST), should be able to detect the void lensing signal with sufficient precision from stacking abun- 
dant medium-sized voids, thus providing direct constraints on the matter density profile of voids independent 
of assumptions on galaxy bias. 

Subject headings: cosmology - large scale structure - dark matter - gravitational lensing: weak 
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1. INTRODUCTION 

As revealed by galaxy redshift surveys, the large-scale 
structure of the universe hierarchically grows into a complex 
network of filaments and galaxy clusters, separated by large, 
nearly empty voids in between. While the bound structures 
contain most of the mass in the Universe, most of the cosmic 
volume is filled with voids. These voids originate from local 
minima of the primordial density field and expand faster than 
the average density of the Universe, clearing out their central 
regions and building up filaments at their boundaries. 

Just as the more studied galaxy clusters, the formation and 
evolution of voids probe the extreme tails of the cosmic den- 
sity distribution, and provide insights into initial conditions 
and the expansion history of the Universe. For example, the 
abundance of massive clusters and large voids increases in 
f{R) gravity models (Li et al. 2012), and non-Gaussian initial 
conditions change the void and cluster abundance in opposite 
directions (Kamionkowski et al. 2009). 

Besides, as the dynamics of underdense regions is domi- 
nated by dark energy at earlier times than the rest of the uni- 
verse, the potential of voids to probe the nature of dark energy 
has been noted recently. Lavaux & Wandelt (201 1) proposed 
an Alcock-Paczynski test on the average shape of stacked 
voids and found that it may outperform the Baryonic Acous- 
tic Oscillation effect as a cosmological probe for a Euclid- 
like survey. In addition, the distribution of void ellipticities 
has been demonstrated to be a powerful probe of the equation 
of state of dark energy (Park & Lee 2007; Lavaux & Wandelt 
2010; Bos et al. 2012). Similarly, the emptiness of voids and 
its evolution over cosmic time can probe the expansion history 
of the universe and modified gravity (e.g., Farrar & Peebles 
2004; Nusser et al. 2005; Peebles & Nusser 2010). 

Observationally, voids are identified in the distribution of 
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galaxies (e.g., Sutter etal. 2012, , and references therein), 
and the matter density profile of voids is inferred from the 
galaxy distribution in and around voids. The latter requires 
assumptions on galaxy biasing, which is not precisely known 
and subject to uncertainties in the efficiency of galaxy forma- 
tion and evolution in underdense environments. 

In this letter we demonstrate that the matter density pro- 
file of voids can instead be constrained directly from the 
average lensing signal of stacked voids. While weak lens- 
ing by voids was previously discussed by Amendola et al. 
(1999), these authors focused on the tangential shear signals 
by individual voids and concluded that void radii larger than 
~ lOOMpc/i"' were required to detect the effect. Compared 
to their analysis, stacking allows us to achieve much higher 
significance in measuring the weak lensing signals induced 
by more abundant, smaller voids, large number of voids will 
be available from future spectroscopic galaxy redshift sur- 
veys, such as the proposed Prime Focus Spectrograph (PES; 
Ellis et al. 2012) instrument on the Subaru telescope and the 
Euclid satellite; while several future imaging surveys, e.g., 
the on-going Subaru Hyper Suprime Cam (HSC), Euclid, and 
the LSST projects, will allow us to measure the stacked weak 
lensing signals by voids. Such direct measurements of void 
density profiles enable cosmological constraints without the 
uncertainties from galaxy biasing models, and may addition- 
ally help in studying the environment dependence of galaxy 
evolution. 



2. METHODS AND MODELS 

We now forecast the sensitivity of stacked weak lensing of 
voids and its ability to measure the average density profile of 
voids. We consider voids in the redshift range of 0.4 < zi < 
0.6, assuming the void catalog is available in a survey area of 
A - 5000 square degrees. We also assume a DETF Stage 
IV deep imaging survey (Albrecht et al. 2006) of the same 
area for the lensing analysis; the mean background density 
of galaxies is assumed to be Wgai = 12 arcmin"-, with a mean 
redshift of Zs ~ 1. Additionally, we consider two specific 
upcoming surveys: the combination of HSC and PES, known 
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Figure 1. Left: void density profile. Right: projected void surface density 
scaled b\R\j. , 

as the SuMIRe survey , with A - 1500 square degrees, Wgai - 
20 arcmin"^, Zs ~ 1, and a Euclid/WFIRST like survey with 
A = 15000 square degrees, Wgai - 30 arcmin"^, ~ 1.2. 

Throughout the paper we use the cosmological parame- 
ters derived from the Wilkinson Microwave Anisotropy Probe 
(WMAP) 7-year results (Komatsu et al. 201 1). 

2.1. Void Models 

While voids are intrinsically ellipsoidal, we assume the av- 
eraged void density profile after stacking to be spherically 
symmetric. In our basic model, the void density profile and 
abundance follow the simulation results of Lavaux & Wandelt 
(201 1). Specifically, we use their fitting formula 



p(r, R^/) = p 



An + A3 1 — 



0.27 -H 0.61 1 — 



3^ 



(1) 



to describe the void density profile within r - Ry, the void 
radius; here p is the mean cosmic matter density, and where 
we have ignored the scatter in the best fit parameters, Aq - 
0.27 ± 0.01 and A3 - 0.61 + 0.03, which has negigable impact 
on our signal-to-noise calculation, assuming that the scatter is 
close to symmetric. 

Since there are few constraints on the outer density profile 
of voids, we choose a continuation outside Ry with a constant 
density wall such that profile is compensated within 2Ry, see 
Fig. 1, in order to model the extended, marginally overdense 
wall structures found to surround voids in recent simulations 
(e.g., Colberg et al. 2008; Lavaux & Wandelt 201 1) and void 
catalogs (Sutter et al. 2012). We explore in Section 4 a range 
of other models, including voids surrounded by steep, over- 
dense ridges, as found in earlier simulations and void catalogs 
(e.g.. Pan et al. 2012). The comoving number density of voids 
in a IMpc/z"' wide size bin is given by (Lavaux & Wandelt 
2011) 

n^oidiRy &[r,r+ IMpc/Zi]) 



(Mpc//z)3 



= 3.5x10"' exp -0.632 



Mpc/Zi 
(2) 



2.2. Weak Lensing Signals 

For a lens at redshift z\ the projected surface mass density S 
is given by 

m Ry) = J dr, {p ( yjieD^izi)? + rl Ry^ - p) , (3) 

where Da{z) is the proper angular diameter distance. In 
the linear regime, which is an excellent approximation for 
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the very weak lensing caused by underdense regions, the 
mean tangential shear yt and magnification fi profile of an az- 
imuthally symmetric mass distribution are given by 

AS(6») 1.(6) 
<rt) (0) = -r^, n(9) ^ 1 + 2k(6) = 1 + 2^ , (4) 



where k is the lensing convergence and AE is the excess sur- 
face mass density defined as 



AS(0) = £(0)-i(0) = 5;erit<rt>(0). 



(5) 



Here Eait is the critical surface mass density weighted by the 
source galaxy distribution dns/dzs. 
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where we have neglected the narrow redshift distribution of 
the lens population and use the mean redshift, zi - 0.5, in 
practice. While Bolejko et al. (2012) recendy pointed out that 
the relativistic Doppler term arrising from correlated pecu- 
liar velocities can have a significant effect on the magnifica- 
tion signal of source galaxies closely associated with the lens- 
ing void, this eff'ect is negligable for our calculation due to 
the large separation between sources and voids, and the large 
number of voids used in the stacking process. 

Observationally, the tangential shear is estimated as the 
mean tangential ellipticity of source galaxies with respect to 
the void center (e.g., Bartelmann & Schneider 2001). A num- 
ber of techniques have been proposed to measure the magni- 
fication profile of galaxy groups and clusters; they broadly 
fall into two categories: magnification bias measurements, 
which utilize the magnification-induced variations in the ap- 
parent number density of background galaxy populations (e.g. 
Hildebrandt et al. 201 1; Umetsu et al. 201 1; Ford et al. 2012; 
Umetsu et al. 2012), and magnification estimators based on 
the change in the apparent sizes and fluxes of individual galax- 
ies (e.g., Huff^ & Graves 201 1; Schmidt et al. 2012). Each of 
these estimators is affected by a variety of systematic effects, 
as described in the references above; a detailed modelling of 
these effects is beyond the scope of this paper, and we will 
absorb this net irreducible observational noise into a simple 
one-parameter model, as described below. 

In order to increase the signal-to-noise ratio (S/N) of the 
lensing observables, the lensing signals are measured by 
stacking large numbers of voids with similar sizes. In prac- 
tice, the size measurement of an individual void in galaxy 
redshift surveys will be very noisy due to sparse sampling in 
the underdense region. Hence, we follow the stacking pro- 
cedure of Lavaux & Wandelt (201 1) and Sutter et al. (2012), 
who stack all voids within a size bin on their barycenters with- 
out rescaling individual voids. Ultimately, the uncertainties 
in the weak-lensing profiles are a combination of observa- 
tional errors and cosmic noise due to projection effects of un- 
corrected large-scale structure along the line of sight. The 
shear (magnification) profile of voids can also be written as 
the small-scale cross correlation function between void cen- 
ters and shear (magnification), and we model their covariance 
matrices in analogy to galaxy-galaxy lensing (e.g. Jeong et al. 
2009). 

The Gaussian covariance of the angular void-shear cross 
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spectrum, C^'''(0 - €^"{1) is given by 



A(2l + 1) 
1 



(c^''(0)' + c^^(c«(o + ^) 



+ — C^''(0 + — , (7) 

with €""{1) the angular cosmic shear power spectrum, cr^ 
the ellipticity dispersion, «v the projected number density of 
voids, and C^^{1) the angular void power spectrum. Note that 
the last term in Eq. 7 is the standard covariance estimate for 
stacked profiles (Hoekstra 2003). In comparison, the corre- 
lation function covariance has additional contributions from 
the void profile {€^"{1)) and from void correlations (C^^(O), 
which increase the statistical uncertainty. The latter includes 
the effective increase in shape noise caused by the overlap of 
projected voids. For the purpose of this S/N estimate we as- 
sume voids to be unclustered; then C^'^ is given by the Fourier 
transform of the void convergence profile, and C^^ is deter- 
mined by its "one-halo" term (see Cooray & Sheth 2002, for 
a review of the halo model). While the void density profiles 
used in the lensing calculations are compensated within 2Ry, 
we will model projected voids as uniform disks with radius 
20V in the calculation of the void clustering power spectrum 
to maximize the potential degradation caused by the projected 
overlap of different voids, and to avoid complications with the 
normalization of the halo model. Then the one-halo term of 
voids with angular radius 6y and angular number density ng^is 
given by 



1 l Jl(219y) 
Wflv I '^v 
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profile is given by 

Cov(M0i,A0),ju(02,A^?)) = 2^ J 



(8) 



with J„ the «-th order Bessel function of the first kind. For 
an ensemble of voids with angular radii in the range dy e 
[0v min, dy max] and radius distribution dnviOy) / dOy it is given 
by' 



(9) 



Finally, we approximate the covariance of the angular shear 
profile in bins [0i - A0/2, 9, + A0/2] with uniform bin-width 
A0 as 



(10) 



where the Fourier transform is evaluated at the area weighted 
bin centers 0,. 

Analogously, the covariance of the angular magnification 
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In the following discussion, we will refer to the two terms in 
Eqs. ( 10, 11 ) proportional to C" as large-scale structure noise, 
and to those proportional to cr^/fj as observational noise. The 
contribution from the void lensing profile {C^") to the covari- 
ance is negligible. 
Due to the high redshift of source galaxies, we assume 
= 0.3 per shear component. The observational noise pa- 
rameter a-fj for the magnification measurement depends on the 
estimator under consideration; for magnification-bias based 

1 /2 

estimators, cr^ - 2l{qf^ ' ), where is the fraction of source 
galaxies used in the magnification measurement, and q ~ l- 
1 .5 is a parameter describing the strength of the magnification 
effect related to the slopes of the galaxy luminosity and size 

distribution functions (Schmidt et al. 2009). For a scaUng re- 

1 y 

lation based estimator (Huff & Graves 201 1), cr^ - 2cr^//s 
with (Tit the scatter in the convergence estimator caused by the 
scatter in the fundamental plane, /s is the fraction of source 
galaxies to which the magnification estimator can be applied, 
and the factor of two propagates from convergence to magnifi- 
cation error. For definiteness, we use cr^ = 2, but note that the 
observational noise in the magnification estimate dominates 
over the large-scale structure noise, so that the error bars are 
proportional to cr^. 

3. RESULTS 

Fig. 2 shows the results of our calculations for the stacked 
void lensing profiles in different size bins, where we have as- 
sumed that the distribution of void sizes within each bin fol- 
lows the form of Eq. (2). This calculation indicates that shear 
and magnification will provide strong constraints on the radial 
matter density profile of small to medium sized voids, if these 
can be identified from a galaxy redshift survey. 

We also include the degradation of the lensing signal by 
stacking voids on incorrectly chosen centers, which is illus- 
trated by the dashed lines in Fig. 2. Specifically, we assume a 
two-dimensional Gaussian distribution of off-center positions 
with variance 1 Mpc/j"', independent of void size. Note that 
this illustrates the effect of severe mis-centering, as it cor- 
responds to a rms offset of 24% of the mean void size for 
the smalllest void size bin ([5, 10] Mpc/h). Due to the shal- 
low, slowly varying density of voids, off-centering has a much 
smaller impact on void lensing profiles than on cluster lensing 
profiles. 

In Table 1 we list the radially integrated S/N for our fidu- 
cial survey, as well as for the SuMIRe and Euclid surveys. 
Here we have included the cosmic covariance between differ- 
ent radial bins, but ignored the covariance between shear and 
magnification measurements caused by large-scale structure 
as the magnification measurements are strongly dominated by 
observational noise. 

The middle and right panels show the signal and noise 
expected for voids in the radius range [5, 10] Mpc//i and 
[10, 15] Mpc/h. As can be seen from Table 1, the integrated 



4 



The Gravitational Lensing Signal of Stacked Voids 



o 

Q. 



0.4 

0.2 

0.0 

-0.2 

-0.4 

-0.6 

-0.8 
-1.0 




0.0 0.5 1.0 1.5 2.0 2.5 

W<6»„> 




0.0 0.5 1.0 1.5 2.0 2.5 

6'/<6»„> 



2.0 2.5 




< 




0.0 0.5 1.0 1.5 2.0 2.5 

6»/<6'„> 

Figure 2. Void lensing profiles obtained from stacking voids in different size bins. The top panel shows the projected surface mass density S, which is proportional 
to the magnification signal. The lower panel shows AZ (Eq. (5)), which is proportional to the mean tangential shear. The inner error bars indicate the cosmic 
noise contiibution due to intervening large scale structure, expected for stacking voids in the redshift range [0.4, 0.6] in size bins as indicated in each plot, and 
the larger error bars also include shape/magnification measurement noise with survey parameters described in Sect. 2. The dashed lines indicate the degradation 
of the lensing signal caused by stacking voids with incorrectly chosen centers, assuming a two-dimensional Gaussian distribution of off-center positions with 
variance lMpc//i. 



S/N decreases for larger voids. In our void model the am- 
plitude of the lensing signal is linearly proportional to void 
size, while the noise increases more rapidly due to the expo- 
nentially decreasing void abundance. Using logarithmic size 
bins for the stacking cannot revert this trend. Hence, for a 
broad range of steeply decreasing void abundance functions, 
the maximum signal-to-noise will be obtained from the small- 
est idendentifiable voids. 

In particular, for void sizes close to the mean spacing be- 
tween galaxies, void finding algorithms may produce incom- 
plete void catalogs. Let c denote the completeness fraction of 
a void catalog, such that c = 1 if all existing voids are identi- 
fied. Then the total S/N scales as 



Table 1 

Void number density and integrated S/N for different void size bins 
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(12) 



with our fiducial {S /N)m from Table 1, which are calculated 
assuming c - 1, and A the survey area. Additionally, the 
S/N for magnification measurements is approximately propor- 
tional to cr„. 



{S/N\ * (S/N) 



(13) 



which holds as long as the error budget is dominated by ob- 
servational noise (cr^, > 0.5). Note that these scaling require 
the redshift distribution of voids and background galaxies to 
stay constant. 



4. DISCUSSION 

The largest theoretical uncertainty of our void model is the 
distribution of matter around voids, that is, the transition from 
underdense void to filament regions. In lack of theoretical 
models for void density profiles we explore how well the lens- 
ing results can differentiate between a range of toy models. 
Specifically, we consider the following alternative models: 

(a) Fitting function of Colberg et al. (2005), continued out- 
side Ry with an extended, low-overdensity ridge. 

(b) Fitting function of Lavaux & Wandelt (2011), contin- 
ued outside Ry, truncated such that the void profile is 
compensated. 
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Figure 3. Impact of void density profile on lensing signal. The left panel shows our fiducial void profile, which is a compensated continuation of the fit function 
of Lavaux & Wandelt (201 1 ) with an extended, low-overdensity ridge, as a thick solid line; the dashed line illustrates a compensated continuation of the same 
central profile with a sharp, high-density ridge; the dash-dotted line shows compensated continuation of the fit function of Colberg et al. (2005) with an extended, 
low-overdensity ridge, and the thin line shows an uncompensated void profile obtained from spherical expansion of an initial top-hat void profile in a ACDM 
universe. The profile illustrated by a dotted line was obtained by extending the ridge of the latter profile by hand. The middle and right panel show the stacked 
magnification and tangential shear signal of voids with different density profiles; en'or bars are the same as described in Fig. 2. 



(c) Spherical top-hat density profile, evolved nu- 
merically in the ACDM cosmology (c.f. 
Sheth & van de Weygaert 2004). 

(d) (c) with an artificially extended rigde. 

Profile (a) is similar to our fiducial model but is less steep, 
having a higher central density. Profiles (b) to (d) have sharp, 
high-density ridges as illustrated in the left panel of Fig. 3. 
The resulting stacked lensing profiles are shown in the middle 
and right panels of Fig. 3. These results illustrate that mag- 
nification and shear signals can measure the extent of wall- 
structures and detect pronounced overdense ridges if these ex- 
ist. However, if individual voids have sharp overdense ridges, 
these ridges will be smoothed out in the stacking process due 
of the ellipsoidal shape of individual voids. Hence to better 
constrain the transition between void and filament regions, in 
practice voids should be stacked along their projected major 
axis instead. 

Additionally, the outer shear profile at r > /?v indicates 
whether the void profile is compensated or underdense: in the 
latter case the the profile falls off like a negative point mass 
outside the angular extent of the void. 

Also note that profiles (b) to (d) have the same central void 
density as our fiducial model. As the void lensing signal is 
affected by the outer void density profile at all projected radii, 
measuring the emptiness of voids from lensing observations 
will require parametric models of the extended void profile. 

5. CONCLUSION 

We have demonstrated that a DETF Stage IV type of galaxy 
redshift survey will allow a clear and robust detection of the 
stacked lensing signal of medium-sized voids with sufficient 
precision, S/A^ > 10, and with S/A^ > 5(15) for the upcom- 
ing SuMIRe (Euclid) survey, when shear and magnification 
are combined, and hence provide strong direct constraints on 
the radial shape of the projected density profile of voids. The 
lensing measurement directly probes the mass distribution in 
and around voids, and provide an unbiased view of these un- 
derdense regions that occupy most of the cosmic volume. It 
will also offer new opportunities to test gravity on cosmic 
scales, but in under-dense regions, thus complementing clus- 
ter based tests (e.g. Lombriseret al. 2012). 
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